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all  three  spatial  coordinates,  and  this  capability  was  demonstrated. 

Two  separate  methods  of  representing  a  sound-speed  profile  (SSP)  based  on  data 
points  were  examined;  Akima  cubic  spline  and  spatial  Fourier  series  (SFS).  The  SFS 
representation  encountered  difficulties  in  accurately  modeling  SSPs.  Various  techniques 
were  applied  to  improve  the  SFS  sound-speed  representation.  While  accurate  sound- 
speed  fits  were  eventually  achieved,  difficulties  remained  in  the  SFS  modeling  of  first  and 
second-order  derivatives  of  the  sound-speed  data. 

The  RRA  algorithm  was  tested  using  the  SFS  sound-speed  representation  and  found 
to  be  significantly  inaccurate. 

A  demonstration  was  conducted  of  the  ability  of  the  SFS  sound-speed  represen¬ 
tation  to  incorporate  randomness  in  the  SSP. 
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I.  INTRODUCTION 


The  primary  focus  of  this  thesis  was  to  test  the  accuracy  and  capabilities  of  a  re¬ 
cursive  ray  acoustics  (RRA)  algorithm  that  can  be  used  to  calculate  the  position,  angles 
of  propagation,  travel  time,  and  path  length  along  a  ray  path  and  to  generate  ray  trace 
plots  for  speeds  of  sound  that  are  functions  of  all  three  spatial  coordinates.  The  RICA 
algorithm  is  an  extension  and  slight  modification  of  an  algorithm  briefly  discussed  by 
Klein  [Ref  1:  pp.  29-31].  The  algorithm  discussed  by  Klein  [Ref  1:  pp.  29-31  ]  was  ex¬ 
pressed  in  terms  of  an  index  of  refraction  while  herein,  the  RR.\  algorithm  is  expressed 
directly  in  terms  of  the  speed  of  sound.  No  computer  simulation  results  were  presented 
by  Klein  [Ref  1:  pp.  29-31  ]. 

The  accuracy  of  the  RRA  algorithm  was  tested  by  comparing  its  computer  simu¬ 
lation  results  with  those  obtained  from  a  previously  developed  ray  acoustics  algorithm 
based  on  solving  a  system  of  first-order  ordinaiy  differential  equations  (ODE)  [Ref  2). 
The  ODE  algorithm  as  applied  to  ray  acoustics  vvas  discussed  in  detail  by  Lim  [Ref  3). 
Both  the  ODE  algorithm  and  the  RRA  algorithm  used  the  IMSL  Version  10  double 
precision  computer  program  DCSAKM  to  obtain  an  Akima  cubic  spline  fit  to  one¬ 
dimensional,  depth-dependent,  sound-speed  data  for  testing  and  comparison  purposes. 
The  computer  simulation  results  are  presented  in  Chapter  III,  Section  A.  An  Akima 
cubic  spline  was  chosen  since  it  is  designed  to  produce  a  curve  that  matches  the  shape 
of  the  data  while  minimizing  oscillations.  In  addition,  the  RRy\  algorithm  was  tested 
and  shown  to  be  capable  of  handling  a  sound-speed  profile  (SSP)  that  is  a  function  of 
all  three  spatial  coordinates.  This  capability  was  demonstrated  by  first  separating  the 
ocean  into  zones  along  the  z  axis.  In  each  zone,  the  SSP  is  a  different  function,  in  gen¬ 
eral,  of  both  cross-range  x  and  depth  y.  A  SSP  with  a  gradient  in  the  cross-range  di¬ 
rection  will  cause  a  sound  ray  to  leave  its  initial  plane  of  propagation.  Chapter  III, 
Section  B,  includes  demonstrations  of  the  special  three-dimensional  capabilities  of  the 
RRA  algorithm. 

A  secondary  objective  of  this  thesis  was  to  exantine  the  use  of  a  spatial  Fourier  series 
(SFS)  to  fit  sound-speed  data  and  to  assess  its  utility  relative  to  an  Akima  cubic  spline 
fit.  The  appeal  of  the  SFS  lay  in  its  representation  of  spectral  content  in  terms  of  spatial 
frequencies,  as  well  as  its  ability  to  represent  derivatives.  The  precision  of  the  two 
sound-speed  fits,  as  well  as  that  of  the  first  and  second-order  derivatives  of  the  sound- 
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speed  data  were  compared  since  the  first-order  derivative  as  well  as  the  sound-speed  fit 
itself  are  utilized  in  the  ray  trace  algoritlims  mentioned  above.  Chapter  IV,  Section  A. 
contains  the  results  of  this  comparison  as  generated  by  computer  simulation.  The  SFS 
representation  of  a  SSP  was  then  employed  within  the  RICA  algorithm  and  tested  by 
using  the  same  SSPs  that  were  used  to  compare  the  ODE  and  RICA  algorithms.  Chapter 
IV.  Section  B,  discusses  the  accuracy  of  the  RICA  algorithm  incorporating  the  SFS  rep¬ 
resentation  of  the  various  SSPs. 

.Additionally,  random  fluctuations  in  the  SSP  were  modeled  using  the  inherent  suit¬ 
ability  of  a  SFS  repre.sentation  for  the  inclusion  of  components  at  particular  harmonics. 
This  random  element  was  intended  to  simulate  fluctuations  generated  by  such  processes 
as  internal  waves  or  by  inaccurate  collection  of  sound-speed  data. 
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II.  THEORETICAL  ANALYSIS 


A.  THE  RRA  ALGORITHM 

If  a  sound  source  is  located  at  =  (-to.j'o.  Zo)  and  if  a  ray  is  launched  with  initial  an¬ 
gles  of  propagation  and  00  (see  Figure  I  ),  then  an  estimate  of  the  position  of  subse¬ 
quent  points  along  the  ray  path  is  given  by  (see  Figure  2  ) 

=  +  Arj.,,  /=  1,2,3,...  (2.1) 

where 

=  ,  (2.2) 

Asi  is  an  increment  of  arc  length,  the  magnitude  of  which  is  at  the  user's  discretion, 

ni  =  UiX  +  Vip  +  Wjz  (2.3) 

is  the  unit  vector  along  the  ray  path  w’here 

Ui  =  sin  Pi  cos  0/ ,  (2.4) 


V,-  =  cos  Pi , 


(2.5) 


and 


Wi  =  sin  Pi  sin  0; 


(2,6) 


are  the  dimensionless  direction  cosines  with  respect  to  the  x,y  and  z  axes,  respectively, 
and  P,  and  0,  are  the  angles  of  propagation. 

Since,  in  general,  is  not  identical  with  from  Klein  [Ref.  1;  pp.  29-31]  we  see 
that,  as  amplified  and  applied  to  the  ocean  acoustics  environment  in  Ziomek  [Ref  4:  pp. 
225  and  230], 


ds 


=  Vrt(r). 


(2.7) 


3 


4 


or,  approximately, 


=  Vn{ri)Asi 


(2.8) 


where  /?(£.)  =  «, , 


«(£/) 


(2.9) 


is  the  dimensionless  index  of  refraction,  c(fi)  is  the  three-dimensional  speed  of  sound  at 

i:q  =  XqX  +yQy  +  ZoZ  ,  (2.10) 

and,  c{r)  is  the  three-dimensional  speed  of  sound  at 

ri  =  xix+y,y  +  z,z  ,  (2.11) 

which  is  the  position  vector  to  a  point  on  the  ray  path  after  the  ith  iteration. 

It  follows  from  Equation  (2.8)  that 


«(£/)«/  -  «(r;_i)«/_i  =  ^ri(ri)Asi 


and,  as  a  result, 


(2.12) 


A  «(r;-i)«/-i  .  ASi  ^  ^ 

- - - +  TT-T  ^"(A-i)  . 


«(a)  «(4) 

where,  for  better  accuracy  in  the  evaluation  of  the  gradient. 


•  .  ^■*<-1  A 

£,-1  =  £j_,  +  , 


(2.13) 


(2.14) 


which  is  the  three-dimensional  spatial  location  half-way  along  the  arc  length  between 
£,_!  and  £,,  and 


n  d  A  5  a  5  a 


(2.15) 


is  the  gradient  expressed  in  rectangular  coordinates.  It  should  be  noted  that  although 
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As,  has  been  shown  as  a  function  of  iteration  number,  it  will  be  constant,  in  general, 
except  when  it  needs  to  be  adjusted  to  properly  handle  surface  and  bottom  reflections 
and  maximum  down  range  limits. 

In  terms  of  the  direction  cosines  u,  v,  and  w. 


«(£,_i)w,_. 

As 

^«(£*_i) 

■  «(/,) 

-  4.  -J±£— 

m) 

CJC 

(2.16) 

7 

7 

s: 

As 

(2.17) 

n{Ld 

■  -h  ■■ 
«(A) 

_  «(r;-i)w/-i  As  Mlj-O 


(2.1S) 


If  c(r,)  =  c(jc„  z,)  is  the  three-dimensional  speed  of  sound,  and  if  Cj  =  cfzo)  is  the 
speed  of  sound  at  the  source  location,  then  it  follows  that 


«(r^-i)  Cq  c(£i)  ^  c(£i) 


(2.19) 


1  Mr) 

n(r,)  ox 


I  c  ii 


dcjr) 

dx  =  C 


(2.20) 


Analogous  reasoning  applied  to  the  y  and  z  coordinates  yields  results  similar  to 
Equation  (2.20) ,  and  substitution  of  these  results  into  Equation  (2.13)  yields 


A  ^(i^l)  A  /  ^  A 


(2.21) 
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In  addition,  although  Klein  did  not  discuss  travel  time  t  and  path  length  s  calcu¬ 
lations,  these  quantities  were  computed  using  the  following  formulae: 


T  = 


(2.22) 


and 


s  =  ^A5,.  (2.23) 

1=1 

Equations  (2.1)  through  (2.23)  constitute  what  is  referred  to  herein  as  the  recursive 
ray  acoustics  (RRA)  algorithm  first  mentioned  in  Chapter  1.  Input  to  this  algorithm  in 
the  form  of  mathematical  formulae  for  the  sound-speed  profile  or  some  particular  nu¬ 
merical  fit  to  sound-speed  data  must  be  provided  to  enable  the  algorithm  to  function. 

B.  SPATIAL  FOURIER  SERIES  REPRESENTATION  OF  SOUND-SPEED  DATA 

Given  an  ocean  D  meters  deep  with  a  one-dimensional,  depth-dependent  sound- 
speed  profile  (SSP),  the  speed  of  sound  c(>')  can  be  represented  by  the  following  spatial 
Fourier  scries  (SFS): 


5 

cO)  =  I  Cq  I  +  2^  I  Cjl  cos{2nfy^y  -  Lc^),  0<y<  D  (2.24) 

j=i 

where  |c,|  and  Lc,  are  the  magnitude  and  phase(angle),  respectively,  of  the  complex 
Fourier  series  coefficient  c,  at  harmonic  s,  |col  is  the  magnitude  of  c,  at  harmonic  s  =  0 
and  is  the  "DC  component"  or  average  value  of  the  sound-speed  data, 

/,,-Tv77  (2-25) 

is  the  fundamental  spatial  frequency  in  the  direction  in  cycles  per  meter  where  N  is  the 
total  number  of  sound-speed  data  points  taken  and  is  the  uniform  sampling  period 
in  the  y  direction  in  meters  per  sample,  and  S  is  the  highest  harmonic.  Derivatives  of 
this  Fourier  series  representation  can  be  directly  obtained. 
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The  complex  Fourier  series  coefficients  c,  are  obtained  by  taking  the  spatial  discrete 
Fourier  transform  (DFT)  of  the  sound-speed  data  as  follows: 


c,  =  ^  Y,<^{nTs)lVs,  s  =  -S . 


(2.26) 


where 


n-\.^cxp{+j2nlX) 


and,  when  sampling  at  the  Nyquist  rate, 


.V=2S-H  1. 


(2.27) 


(2.28) 


An  attempt  to  represent  a  function  by  a  truncated  Fourier  series  will  always  result 
in  oscillations  about  the  functions'  true  value.  These  oscillations  normally  have  small 
enough  amplitudes  that  they  are  of  little  consequence.  However,  an  attempt  to  use  a 
SFS  representation  of  a  linear  SSP  with  a  single  constant  gradient  results  in  significant 
discrepancies  in  amplitude  near  the  ocean  surface  and  bottom,  due  to  the  inherent  dis¬ 
continuity  in  the  SSP  that  exists  at  these  boundaries.  This  difficulty,  known  as  "Gibbs' 
phenomenon",  does  not  diminish  as  more  harmonics,  that  is,  sound-speed  data  points, 
are  added  to  the  Fourier  series. 

C.  Lanezos  [Ref  5]  suggested  that  this  phenomenon  would  be  ameliorated  if  the 
Fourier  series  could  be  forced  to  converge  more  rapidly  by  modifying  the  Fourier  series 
coelficients  c,  by  a  correction  factor,  such  that  the  Lanezos  smoothed  Fourier  series  co¬ 
efficients  are  given  by 


sin(57r/S) 


Cs,  s¥^0, 


(2.29) 


where  c,  is  given  by  Equation  (2.25).  Hamming  [Ref  6]  states  that  S  in  Equation  (2.29) 
is  usually  replaced  with  5  -t-  1  and,  thus,  our  representation  of  the  correction  factor  will 
be 


sin(.srt/(.S  +  1)) 
5rr/(S  +  1) 


q,  5=5^0. 


(2.30) 
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C.  LNCORPORATION  OF  RANDOMNESS  INTO  THE  SOUND-SPEED 
PROFILE 

The  speed  of  sound  as  a  function  of  depth,  f(>)  ,  can  be  represented  as  the  sum  of 
deterministic  and  random  components  as  follows: 

cO-O  =  CnCl’)  +  (2.31) 

with  c^(>)  being  the  deterministic  and  C;,(>)  the  random  speeds  of  sound  as  functions  of 
depth  y. 

The  autocorrelation  function  of  C;,0-)  is  ,  for  any  particulars  and  s'  , 

=  (2.32) 


where  is  the  complex  conjugate  of  Cr  at  depths',  0  <s'  ^  D  . 
If 


.V 

V 


CftCv)  =  2^  exp(  -  j27znfy,y)< 

n=^N' 


for  some  set  of  Fourier  series  cocfTicients  c„  and 


(2.33) 


(2.34) 


for  D  the  ocean  depth  in  meters,  then 

1 

c„  =  -^  J  c,(s)  exp(  +j2nnfy^y)dy.  (2.35) 

If  Cr(s)  is  a  random  function  of  depth,  then  c„  will  be  a  random  variable. 

Since 


=  -y  f  f  F[c/j(s-)4(>-')}  exp(  -1-  J2r:fyj^ny  -  n’y’))dydy'  ,  (2.36) 

D  Jq  •’o 


wiiere 


=  /^crCv./). 


(2.37) 
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and  (c^)'  is  the  complex  conjugate  of  ,  and  if  it  is  assumed  that  Cr(v)  is  a  wide-sense 
stationan’  (  i.e.,  statistically  homogeneous )  random  process,  then 

R,p,y)  =  Rc,m,  (2.38) 


where  A>’  =  y  —  /  •  If  R^f^  is  represented  by  a  finite  Fourier  series 

,v 

=  Yj  “  72?r^M^Ay) ,  (2.39) 

m=-S 

« 

then,  as  developed  in  Thomas  [Ref  7:  pp.  148-153], 

~  2.V-I-  1  (2-40) 

where  is  the  power  spectral  density  of  c^fy). 

Thomas  [Ref  7:  pp  148-153]  goes  on  to  show  that 

E[CnCn]  =  r^^nn-^  (2-41) 

where  <3„.  is  the  Kronecker  delta  function.  Therefore, 

Hence,  if  c„  is  a  zero  mean  Gaussian  random  variable  with  variance  r, ,  that  is,  if 

c„  =  AX0,r,),  (2.43) 

then  Cj,(y)  will  be  a  zero  mean  Gaussian  random  process. 

Therefore, 


=  ^CR(Ay)  =  exp(  -  J2nfy^/S.y),  (2.44) 

n 

which  represents  the  Inverse  Discrete  Fourier  Transform  of  the  power  spectral  density. 

Thus,  knowledge  of  the  power  spectral  density  of  c„(y)  enables  calculation  of  r„  via 
Equation  (2.42),  which  in  turn  yields  c„  via  (2.43)  and  C/,(y)  via  (2.33). 
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III.  COMPUTER  SIMULATION  RESULTS 


A.  RAY  ACOUSTICS  -  COMPARATIVE  ANALYSIS  OF  ODE  VERSUS  RRA 
ALGORITHMS 

Four  distinct  sound-speed  profiles  (SSPs)  were  examined  for  purposes  of  compar¬ 
ative  analysis.  For  all  cases,  a  sound  source  spatial  location  of  jtq  =  0  m  and  z,  =  0  m 
was  used.  The  source  depth  was  established  at  a  point  half-way  between  the  ocean 
surface  and  bottom,  100  m  in  most  cases.  A  constant  ocean  depth  of  200  m  was  mod¬ 
eled,  unless  otherwise  noted,  and  for  all  cases  three  separate  rays  were  examined.  All 
rays  were  propagated  initially  in  the  direction  of  the  positive  z  a.xis,  that  is,  =  90“,  but 
/Jo  was  altered  in  order  to  generate  unique  rays.  For  SSPs  1  to  3,  the  launch  angle  /?, 
was  set  to  45°,  85°,  and  135°  .  SSP  4  was  designed  to  simulate  a  SOFAR  channel. 
Launch  angles  were  therefore  chosen  to  result  in  the  first  two  rays  propagating  down 
range  without  reflection,  i.e.,  =  85°  and  95°,  while  the  third  launch  angle  was  set  to 

=  135°  to  test  reflection  as  well. 

The  four  SSPs  considered  are  defined  as  follows: 

•  SSP  1:  Constant  speed  of  sound  (  c  =  1500  m  s  ) 

•  SSP  2:  Linear  SSP  with  a  positive  gradient  (  c(j)  =  c(0)  -H  gy  m's  where  the  speed 

of  sound  at  the  surface  c(0)  =  c{y  =  0)  =  1.500  m/s  and  g  ==  -f0.017  1,  s ) 

•  SSP  3:  Linear  SSP  with  a  negative  gradient  ( cC^)  =  +  sy  ni  s  w'here  the 

speed  of  sound  at  the  surface  c(0)  =  c(y  =  0)  =  1500  m/s  and  g  =  —0.017  1/s  ) 

•  SSP  4:  Parabolic  SSP  with  vertex  at  >’o  where  c(.Vo)  =  1490  m/s  and  the  speed  of 
sound  at  both  the  ocean  surface  and  bottom  is  equal  to  1500  m  s 

( (■(^■)  =  c(>\,)  -t-  a(y  —  m/s  with  a  =  0.001  for  an  ocean  depth  of  200 

m). 

The  above  four  SSPs  generated  nine  sound-speed  data  points  which  were  input  to 
the  Akima  cubic  spline  routine.  The  results  from  the  ODE  algorithm  were  used 
throughout  as  the  benchmark  for  comparison  purposes.  The  ODE  algorithm  used  the 
IMSL  Version  10  double  precision  computer  program  DIVPBS  in  order  to  solve  the 
system  of  first-order  ordinary  differential  equations.  A  range  step  size  of  2.0  m  and  an 
error  tolerance  of  1  x  lO  ’  was  used  in  DIVPBS,  except  for  SSP  3,  Ray  I  in  Table  3, 
where  the  tolerance  had  to  be  reduced  to  1  x  10  ‘  in  order  for  the  ODE  algorithm  to  run. 
The  arc  length  step  size  used  in  the  RRA  algorithm  was  2.0  m  unless  otherwise  noted. 
The  results  of  the  various  simulation  runs  are  presented  in  the  following  seven  tables. 
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Table  1.  RAY  ACOUSTICS  RESULTS  FOR  SSP  1,  RAYS  I  THROUGH  3 


Method 

J'rj  (m) 

Po  (Jcg) 

r  (km) 

y  (m) 

P  (deg) 

T  (sec) 

s  (km) 

ODE 

100.0 

45.0 

10.0 

100.00 

45.000 

9.428090 

14.1421 

RR/\ 

100.0 

45.0 

10.0 

99.95 

45.000 

9.428091 

14.1421 

ODE 

100.0 

85.0 

10.0 

174.89 

85.000 

6.692132 

10.0382 

RR.\ 

100.0 

85.0 

10.0 

174.88 

85.000 

6.692132 

10.0382 

ODE 

100.0 

135.0 

10.0 

100.00 

135.000 

9.428090 

14.1421 

RRA 

100.0 

135.0 

10.0 

100.05 

135.000 

9.428091 

14.1421 

For  SSP  1  (constant  speed  of  sound),  Table  1  shows  that  at  a  horizontal  range 
r  =  10  km,  the  differences  between  the  two  algorithms  are  as  follows: 

•  between  0.01  and  0.05  m  in  depth  >'  and 

•  I  u  sec  in  travel  time  t. 
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Table  2.  RAY  ACOUSTICS  RESULTS  FOR  SSP  2,  RAYS  1  THROUGH  3 


.Method 

>0  (m) 

(deg) 

r  (km) 

y  (m) 

/?  (deg) 

T  (sec) 

s  (km) 

ODE 

100.0 

45.0 

10.0 

99.99 

45.000 

9.417413 

14.1421 

RR.A 

100.0 

45.0 

10.0 

99.95 

45.000 

9.417422 

14.1421 

ODE 

100.0 

85.0 

10.0 

162.50 

85.487 

6.683492 

10.0374 

RR.\ 

100.0 

85.0 

10.0 

162.57 

85.487 

6.683496 

10.0374 

ODE 

100.0 

135.0 

10.0 

100.01 

135.000 

9.417413 

14.1421 

RR/\ 

100.0 

135.0 

10.0 

100.05 

135.000 

9.417422 

14.1421 

For  SSP  2  (linear  SSP  with  a  positive  gradient),  Table  2  shows  that  at  a  horizontal 
range  r  =  10  km,  the  difierences  between  the  two  algorithms  are  as  follows; 

•  between  0.04  and  0.07  m  in  depth  >•  and 

•  between  4  and  9  ^  sec  in  travel  time  t. 
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Table  3.  RAY  ACOUSTICS  RESULTS  FOR  SSP  3,  RAYS  1  THROUGH  3 


Method 

>b  (m) 

A  (lieg) 

r  (km) 

y  (m) 

ji  (deg) 

T  (sec) 

5  (km) 

ODE 

100.0 

45.0 

10.0 

99.99 

45.000 

9.438784 

14.1421 

RRA 

100.0 

45.0 

lO.O 

99.95 

45.000 

9.438795 

14.1421 

ODE 

100.0 

85.0 

10.0 

168.89 

84.512 

6.699404 

10.0380 

RRA 

100.0 

85.0 

10.0 

168.96 

84.511 

6.699409 

10.0380 

ODE 

100,0 

135.0 

10.0 

100.01 

135.000 

9.438784 

14.1421 

RRA 

100.0 

135.0 

10.0 

100.05 

135.000 

9.438794 

14.1421 

For  SSP  3  (linear  SSP  with  a  negative  gradient),  Table  3  shows  that  at  a  horizontal 
range  r  =  10  km,  the  dilTerences  between  the  two  algorithms  are  as  follows: 

•  between  0.04  and  0.07  m  in  depths  , 

•  0.001  deg  in  angle  of  arrival  /? ,  and 

•  between  5  and  11  m  sec  in  travel  time  r. 

The  data  presented  in  Tables  1  through  3  demonstrate  that  the  two  algorithms  are 
in  very-  good  agreement.  Note  that  the  accuracy  of  the  RRA  algorithm  can  always  be 
improved  by  reducing  the  arc  length  step  size.  This  is  demonstrated  in  Tables  5  and  6. 


14 


Table  4.  RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  I  THROUGH  3 


Method 

To  (m) 

/^o  (Jeg) 

r  (km) 

y  (m) 

/?  (deg) 

z  (sec) 

s  (km) 

ODE 

100.0 

85.0 

10.0 

37.06 

87.244 

6.710387 

10.0183 

RIGA 

100.0 

85.0 

10.0 

36.63 

87.224 

6.710335 

10.0184 

ODE 

100.0 

95.0 

10.0 

162.94 

92.756 

6.710387 

10.0183 

RR.\ 

100.0 

95.0 

10.0 

163.37 

92.776 

6.710335 

10.0184 

ODE 

100.0 

135.0 

10.0 

144.9" 

134.922 

9.448852 

14.1104 

RR.\ 

_ 

100.0 

135.0 

lO.O 

144.24 

134.930 

9.449238 

14.1110 

For  SSP  4  (parabolic  SSP),  Table  4  shows  that  at  a  horizontal  range  r  =  10  km,  with 
the  ocean  bottom  at  200  m,  the  difTerences  between  the  two  algorithms  are  as  follows: 

•  between  0.43  and  0.73  m  in  depth  y  , 

•  between  0.008  and  0.02  deg  in  angle  of  arrh'al  ^  , 

•  between  52  and  386  ix  sec  in  travel  time  t  ,  and 

•  between  0.1  and  0.6  m  in  path  length  5. 

As  is  indicated  in  Table  4,  the  differences  in  travel  time  are  relatively  large  when 
compared  with  previous  cases.  This  lower  level  of  accuracy  was  considered  unacceptable 
since,  in  general,  accuracy  in  the  millisecond  region  is  normally  required  for  travel  time. 
An  attempt  to  improve  accuracy  was  made  by  reducing  the  arc  length  increment  used 
in  the  RRA  algorithm  from  2.0  to  1.0  m  (see  Tables  5  and  6). 

The  parabolic  SSP  which  produced  Table  4  has  the  most  extreme  gradient  of  all  the 
SSPs  examined  and,  thus,  the  two  rays  launched  at  85°  and  95°  proceed  down  range 
without  reflection  from  the  surface  or  bottom,  as  illustrated  in  Figures  3  and  4.  This 
extreme  gradient  was  considered  to  have  been  the  probable  cause  of  the  discrepancies 
noted  above  and,  therefore,  the  gradient  was  reduced  and  the  ocean  depth  increased  to 
permit  an  examination  of  a  similar  profile  with  a  more  typical  gradient  (see  Table  7). 
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Table  5.  RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  2.0 
M  ARC  LENGTH  STEP  SIZE 


Method 

Jo  (m) 

Po  (deg) 

r  (km) 

j(m) 

P  (deg) 

T  (sec) 

s  (km) 

ODE 

100.0 

85.0 

5.0 

64.25 

85.597 

3  354779 

5.0088 

RRA 

100.0 

85.0 

5.0 

64.13 

85.581 

3.354757 

5.0088 

ODE 

100.0 

95.0 

5.0 

135.75 

94,403 

3.354779 

5.0088 

RRA 

100.0 

95.0 

5.0 

135.87 

94.419 

3.354757 

5.0088 

ODE 

100.0 

IJ5.0 

5.0 

77.50 

45.019 

4.724411 

7.0552 

RRA 

100.0 

135.0 

5.0 

11.66 

45.017 

4.724506 

7.0553 

The  data  in  Table  5  for  a  horizontal  range  of  5  km  and  using  an  arc  length  step  size 
of  2.0  m  was  generated  in  order  to  enable  a  direct  comparison  with  the  arc  length  step 
size  in  the  RRA  algorithm  set  to  1.0  m,  as  shown  in  Table  6.  The  reduction  in  final 
range  from  10  to  5  km  was  necessitated  by  the  structure  of  the  RRA  computer  code, 
which  required  large  quantities  of  computer  storage  space  to  be  allocated,  particularly 
as  arc  length  step  size  was  decreased.  This  drawback  could  be  rectified  by  straightfor¬ 
ward  modifications  to  the  simulation  program. 

For  SSP  4  (parabolic  SSP),  Table  5  shows  that  at  a  horizontal  range  of  r  =  5  km, 
the  differences  between  the  two  algorithms  are  as  follows: 

•  between  0.12  and  0.16  m  in  depths' , 

•  between  0.002  and  0.016  deg  in  angle  of  arrival  P  , 

•  between  22  and  95  n  sec  in  travel  time  t  ,  and 

•  0. 1  m  in  path  length  s 
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Table  6.  RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  1.0 
M  ARC  LENGTH  STEP  SIZE 


Method 

Vo  (m  ) 

Po  (deg) 

r  (km) 

y  (ni) . 

P  (deg) 

T  (sec) 

s  (km) 

ODE 

100.0 

85.0 

5.0 

64.25 

85.597 

3.354779 

5.0088 

RICA 

100.0 

85.0 

5.0 

64.19 

85.589 

3.354768 

5.0088 

ODE 

100.0 

95.0 

5.0 

135.75 

94.403 

3.354779 

5.0088 

RR.A 

100.0 

95.0 

5.0 

135.81 

94.411 

3.354768 

5.0088 

ODE 

100.0 

135.0 

5.0 

77.50 

45.019 

4.724411 

7.0552 

RRA 

100.0 

135.0 

5.0 

77.56 

45.018 

4.724458 

7.0553 

For  SSP  4  (parabolic  SSP),  Table  6  shows  that  at  a  horizontal  range  of  r  =  5  km, 
and  with  the  RR.A  arc  length  step  size  set  to  1.0  m,  the  dilTercnces  between  the  two  al¬ 
gorithms  are  as  follows: 

•  0.06  m  in  depth  , 

•  between  0.001  and  0.008  deg  in  angle  of  arrival  p  , 

•  between  1 1  and  47  n  sec  in  travel  time  t  ,  and 

•  0.1  m  in  path  length  s. 

A  significant  improvement  in  the  consistency  of  travel  time  calculations  can  be  seen 
by  comparing  Tables  5  and  6.  Cutting  the  arc  length  increment  in  half  reduced  the  travel 
time  discrepancy  by  approximately  one-half.  Further  reductions  in  arc  length  increment, 
an  input  parameter  to  the  simulation  program,  w'ould  yield  further  improvements  in 
travel  time  accuracy. 
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Table  7.  RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  2.0 
M  ARC  LENGTH  STEP  SIZE,  BOTTOM  DEPTH  1000  M 


Method 

>'o  (m) 

Po  (deg) 

r  (km) 

>'(m) 

A (deg) 

T  (sec) 

s  (km) 

ODE 

500.0 

85.0 

10.0 

777.41 

93.388 

6.705898 

10.0150 

RRA 

500.0 

85.0 

10.0 

177.48 

93.389 

6.705893 

10.0150 

ODE 

500.0 

95.0 

10.0 

222.59 

86.612 

6.705898 

10.0150 

RRA 

500.0 

95.0 

10.0 

222.52 

86.611 

6.705893 

10.0150 

ODE 

500.0 

135.0 

10.0 

545.01 

134.997 

9.448814 

14.1104 

RRA 

_ 

500.0 

135.0 

10.0 

545.00 

134.997 

9.448829 

14.1104 

For  SSP  4  (parabolic  SSP)  with  vertex  c(500)  =  1490  m/s,  speed  of  sound  at  the 
ocean  surface  c(0)  ==  1500  m/s,  and  speed  of  sound  at  the  ocean  bottom  c(lOOO)  =  1500 
m/s,  i.e.,  an  ocean  depth  of  1000  m;  with  c{y)  =  c(500)  +  a(y  -  500^  where 
a  =  0.00004,  at  a  horizontal  range  r  =  10  km,  the  differences  between  the  two  algo¬ 
rithms  are  as  follows: 

•  between  0.01  and  0.07  m  in  depths, 

•  0.001  deg  in  angle  of  arrival  ^  and, 

•  between  5  and  15  ^  sec  in  travel  time  t. 

Table  7,  when  compared  with  Table  4,  both  pertaining  to  an  arc  length  step  size  of 
2.0  m,  demonstrates  that  a  more  realistic  gradient  of  0.02  1/sec  enables  the  RRA  algo¬ 
rithm  to  handle  propagation  in  a  SOFAR  channel  to  a  range  of  10  km  with  a  much 
improved  degree  of  accuracy.  Of  course,  this  could  be  improved  still  further  by  reducing 
the  arc  length  step  size. 


20 


B.  RRA  THREE-DIMENSIONAL  CAPABILITIES 

Figure  5  demonstrates  the  ability  of  the  RICA  algorithm  to  handle  launch  angles 
other  than  =  90°  .  This  figure  was  produced  using  the  parabolic  SSP  4  with 

cO)  =  c(>'o)  -I-  (7fv  -  Vo)',  0  <  >-<  200m,  (3.1) 

where  Jo  =  100  m,  c(io)  =  1490.0  m,'s,  a  =  0.001  -p—  ,  and  0o=  85°  .  Note  that  the 
ray  stays  in  the  plane  created  by  this  angle,  as  illustrated  by  the  projection  of  the  ray 
onto  the  A'  Z  plane  as  a  straight  line. 

The  RRA  algorithm  also  has  the  capability  to  divide  the  ocean  into  separate  zones 
as  a  function  of  z.  Each  zone  can  have  its  own  SSP  as  a  function  of  both  cross  range  x 
and  depth  j;.  For  computer  simulation  purposes,  the  speed  of  sound  was  defined  as  fol¬ 
lows; 

c{x,  y)  =  cCvo)  -H  a(y  -  j’o)^  -I-  Jf ,  0  <  y  <  200  m  , 

where  yo  =  100  m,  c(yo)  =  1490.0  m  s,  a  ~  0.001  and 

0.0  l,s,  2  5  km 

-0.05  1/s,  5  <  2  ^  10  km 

The  particular  gradient  in  the  x  direction  was  selected  solely  for  computer  simulation 
purposes  to  illustrate  the  capability  of  the  RRA  algorithm  to  accommodate  sound  speed 
as  a  function  of  both  ar  and  y.  Figure  6  illustrates  the  ray  trace  provided  by  the  RRA 
algorithm  for  the  SSP  given  by  Equations  3.2  and  3.3,  showing  the  ray  staying  in  the 
Y  Z  plane  until  reaching  5  km,  where  the  gradient  in  the  x  direction  comes  into  efi'ect, 
resulting  in  out-of-plane  propagation. 


(3.2) 


(3.3) 
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RECURSIVE  RflY  RCOUSTICS 


22 


RECURSIVE  RflY  ACOUSTICS 


Figure  6.  Ray  Trace  Based  on  a  SSP  as  a  Function  of  x,  y  and  z  > 
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IV.  REPRESENTATION  OF  SOUND-SPEED  DATA 


A.  COMPARATIVE  ANALYSIS  -  SPATIAL  FOURIER  SERIES  VERSUS  AKIMA 
CUBIC  SPLINE  REPRESENTATIONS  OF  SOUND-SPEED  DATA 

Five  dilTerent  exact  mathematical  functions  were  used  to  generate  1 1  sound-speed 
data  points  evenly  spaced  in  depth  in  a  200  m  deep  ocean.  Akima  cubic  spline  and 
spatial  Fourier  series  representations  of  the  SSPs  were  produced  from  these  data  points. 
Percent  error  calculations  of  their  sound-speed  fits  as  well  as  their  first  and  second-order 
derivative  fits  were  performed  with  respect  to  the  known  exact  values.  If  the  error  be¬ 
tween  the  exact  and  modeled  values  was  less  than  10"‘ ,  percent  error  was  set  to  zero. 

The  exact  functions  examined  were  as  follows: 

•  Constant  speed  of  sound  (  c  =  1500  m.'s  ) 

•  Linear  SSP  with  a  positive  gradient  (  c(>’)  =  1500.0  -I-  O.OlTj;  m/s  ) 

•  Linear  SSP  with  a  negative  gradient  (  c(>')  =  1500.0  —  0.017_y  m/s  ) 

•  Half-period  sine  wave  SSP  (  c(v)  =  1500.0  —  25.0  sin(  )  m/s  ) 

•  Parabolic  SSP  (  c(>)  =  1-450.0  -I-  0.005(y  —  lOO.Oy  m/s ) 

Tables  8  and  9  with  Figures  7,  8  and  9  illustrate  the  difficulties  that  the  Fourier 
representation  encounters  in  dealing  with  data  derived  from  an  exact  constant  speed  of 
sound.  No  noticeable  difference  in  the  sound-speed  ,.lt  itself  is  observed  for  either  the 
Akima  cubic  spline  or  the  Fourier  series  representation,  but  the  latter  generates  some 
small  errors  in  both  the  first  and  second-order  derivatives. 

Tables  10  and  11  along  with  Figures  10,  11  and  12  compare  the  Akima  cubic  spline 
and  Fourier  series  representations  of  sound-speed  data  generated  from  a  linear  SSP  with 
a  positive  gradient.  In  Figure  10  and  Table  II  it  can  be  seen  that  the  Fourier  represen¬ 
tation  does  not  precisely  match  the  sound-speed  function,  although  Figure  10  and  Table 
10  show  that  the  Akima  cubic  spline  does.  Problems  similar  to  those  noted  for  a  con¬ 
stant  speed  of  sound  can  be  observed  by  examining  Figures  11  and  12  and  Table  11. 
Discrepancies  in  the  derivatives  are  now  much  larger. 

Tables  12  and  13  along  with  Figures  13,  14  and  15  show  analogous  results  for  a 
linear  SSP  with  a  negative  gradient.  Once  again,  neither  the  sound-speed  fit  nor  the 
derivatives  are  precisely  modeled  by  the  Fourier  series  representation. 
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Table  8.  AKIMA  CUBIC  SPLINE  RESULTS  -  CON¬ 
STANT  SPEED  OF  SOUND 
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Table  9.  FOURIER  SERIES 
SPEED  OF  SOUND 


RESULTS  -  CONSTANT 
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TECHNIQUES  FOR  NUMERICAL  INTERPOLATION  OF  C(Y) 


TECHNIQUES  FOR  NUMERICnt  INTERPOLATION  OF  C(Y) 
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TECHNIQUES  FOR  NUMERICAL  INTERPOLATION  OF  C(Y) 


Table  10.  AKIMA  CUBIC  SPLINE  RESULTS  -  LINEAR 
SSP  VVn  H  A  POSITIVE  GRADIENT 
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LINEAR  SOUND-SPEED  PROFILE  HITH  A  POSITIVE  GRADIENT 


Table  II.  FOURIER  SERIES  RESULTS  -  LINEAR  SSP 
WITH  A  POSITIVE  ORADIENT 
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Figure  11.  First-Order  Derivative  Fit  -  Linear  SSP  with  a  Positive  Gradient. 
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LINEAR  SOUND-SPEED  PROFILE  KITH  A  NEGATIVE  GRADIENT 


Table  12.  AKIMA  CUBIC  SPLINE  RESULTS  -  LINEAR 
SSP  WITH  A  NEGATIVE  GRADIENT 
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Table  13.  FOURIER  SERIES  RESULTS  -  LINEAR  SSP 
WITH  A  NEGATIVE  GRADIENT 
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Figure  14.  First-Order  Derivative  Fit  -  Linear  SSP  with  a  Negative  Gradient. 
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MERICflL  INTERPOLATION  OF  C(Y) 


Figure  15.  Second-Order  Derivative  Fit  -  Linear  SSP  with  a  Negative  Gradient. 
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The  above  three  exact  sound-speed  functions  were  used  to  compare  the  Fourier  se¬ 
ries  representation  with  the  Akima  cubic  spline  since  the  Fourier  fit  can  be  expected  to 
be  at  its  worst  since  we  are  try  ing  to  form  linear  functions  from  trigonometric  ones.  The 
next  two  functions  examined,  however,  are  non-linear;  the  half-period  sine  wave  and  the 
parabola. 

Tables  14  and  15  along  with  Figures  16,  17  and  18  compare  the  Akima  cubic  spline 
and  Fourier  series  representations  of  sound-speed  data  generated  from  a  half-period  sine 
wave.  Both  techniques  encounter  some  dilFiculty  here,  with  the  Fourier  representation 
generally  behaving  better  near  the  midpoint  of  the  data  than  the  Akima  cubic  spline  fit 
for  sound  speed  and  both  derivatives.  However,  near  the  surface  and  bottom,  the  errors 
in  the  Fourier  representation  become  much  larger  than  those  of  the  Akima  cubic  spline 
fit. 

Tables  16  and  17  along  with  Figures  19,  20  and  21  eompare  the  Akima  cubic  spline 
and  Fourier  scries  representations  of  sound-speed  data  generated  from  a  parabola. 
Surprisingly,  the  Akima  cubic  spline  representation  matches  sound  speed  and  derivatives 
very  well.  The  Fourier  representation  behaves  in  a  manner  siinilar  to  its  handling  of  the 
half-  period  sine  wave.  It  produces  small  errors  near  the  midpoint  of  the  data  and  large 
discrepancies  at  either  end  for  sound  speed  as  well  as  in  the  first  and  second-order  de¬ 
rivatives. 

The  disappointing  performance  of  the  Fourier  series  representation  of  sound  speed 
and  of  the  first  and  second-order  derivatives  of  sound  speed  led  to  various  attempts  to 
salvage  it.  First,  Lanezos  smoothing  as  discussed  in  Chapter  II.  Section  C,  was  imple¬ 
mented.  This  ameliorated  the  oscillations  of  the  Fourier  representation  as  shown  in 
Tables  18  and  19  and  Figures  22  through  27  for  the  half-period  sine  wave  and  parabolic 
SSPs.  It  can  be  seen  that  Lanezos  smoothing  results  in  a  smoother  fit  for  sound  speed 
and  for  the  derivatives  to  a  point  closer  to  the  end  points  of  the  data  but,  nevertheless, 
leaves  significant  discrepancies  near  the  end  points. 
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Table  14.  AKIMA  CUBIC  SPLINE  RESULTS 
HALF-PERIOD  SINE  WAVE  SSP 
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Table  15.  FOURIER  SERIES  RESULTS  -  HALF-PERIOD 
SINE  WAVE  SSP 
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TECHNIQUES  FOR  NUHERICRI.  INTERPOLATION  OF  C(Y) 


WITHOUT  UtCZOS  SnOCTHlNS 


TECHNIQUES  FOR  NUMERICAL  INTERPOLATION  OF  C(T) 


Figure  17.  First-Order  Derivative  Fit  -  ({alf-Pcriod  Sine  Wave  SSP. 
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Table  16.  AKIMA  CUBIC 
PARABOLIC  SSP 


SPLINE  RESULTS 
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Table  17.  FOURIER  SERIES  RESULTS  -  PARABOLIC 
SSP 
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DUES  FOR  NUMERICfiL  INTERPOLATION  OF  C(Y) 


ORieiNfL  NUMBER  or  Dnn  POINTS!  11  NUraCR  or  SYNTICTIC  ORTA  POINTS  BETORC  RM)  nriER  1 0 


TECHNIQUCS  FOR  NUHERICRL  INTERPOLATION  OF  CIY) 


Figure  20.  First-Order  Derivative  Fit  -  Parabolic  SSP. 
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ORiaiHn.  NursER  or  ORin  roiNni  ii  nuhbcr  or  snmciic  omn  points  betorc  iwd  RnEitio 

WITHOUT  LAOS  SHOOTHINa 


Table  18.  FOURIER  SERIES  RESULTS  WITH  LANCZOS 
SMOOTHING  -  HALF-PERIOD  SINE  WAVE 
SSP 
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Table  19.  FOURIER  SERIES  RESULTS  WITH  LANCZOS 
_  SMOOTHING  -  PARABOLIC  SSP  . 
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ECHNIQUES  FOR  NUMERICAL  INTERPOLPTION  OF  CtY) 


Figure  22.  Sound  Speed  Fit  with  Lanczos  Smoothing  -  Half- Period  Sine  Wave 
SSP. 
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cniHioous  sazoKn  kiih 


TECHNIQUES  FOR  NUMERICAL  INTERPOLATION  OF  C(Y) 


Figure  23.  First-Order  Derivative  Fit  with  Lanczos  Smoothing  -  Half-Period  Sine 
Wave  SSP. 
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ORioiNH.  fueot  or  oniH  PoiNTSiU  wteER  or  sYMTHCTic  onin  POINTS  BcroRC  no  nrrcRio 

WITH  UWC2DS  SKlOTHne 


TECHNIQUES  FOR  NUMERICni  INTERPOLRTION  OF  C(Y) 


Figure  24.  Second-Order  Derivative  Fit  with  Lanczos  Smoothing  -  Half-Period  Sine 
Wave  SSP. 
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Him  UWCZOS  SMOOTHINB 


TECHNIQUES  FOR  NUMERICRL  INTERPOLRTION  OF  C(Y) 


Figure  25.  Sound  Speed  Fit  with  Lanczos  Smoothing  -  Parabolic  SSP. 
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WITH  UMCZDS  SmtHtW 


TECHNIQUES  FOR  NUhCRICflL  INTERPDLRTION  OF  C(Y) 


Figure  26.  First-Order  Derivative  Fit  with  Lanczos  Smoothing  -  Parabolic  SSP. 
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TECHNIQUES  FOR  NUMERICRL  INTERPOLRTION  OF  CtY) 


WITH  LRICZOS  skoothim: 


The  difTiculty  in  achieving  a  good  fit  at  or  near  the  end  points  of  the  data  using  a 
spatial  Fourier  series  representation  is  a  result  of  the  instantaneous  commencement  and 
cessation  of  sound-speed  data  at  the  ocean  surface  and  bottom.  Therefore,  an  attempt 
to  overcome  this  dilTiculty  was  made  by  adding  synthetic  data  points  above  the  ocean 
surface  and  below  the  ocean  bottom  so  as  to  enable  the  Fourier  series  to  "settle  down" 
before  the  real  data  begins. 

Synthetic  data  was  created  by  first  examining  the  three  data  points  nearest  the  ocean 
surface  and  bottom.  Similar  procedures  were  used  to  generate  synthetic  data  both  above 
the  su.face  and  below  the  bottom.  Therefore,  only  the  surface  case  will  be  examined. 

Let 


"'10  = 


c(v,)  —  c(0) 


Ti 


where  c(0)  =  ciy  =  0)  ,  and 


'"21  =• 


yi  -  y\ 


where  >■,  is  depth  in  meters  at  the  ith  data  point.  If  we  let 

Am  =  ^2,  -  .-Tio 


and  if 


I  Awl  =  l'”2i  -  '"lol 


for  some  arbitrarily  small  e,  then  linear  extrapolation  is  used,  i.e., 
c(j-;)  =  c(0)  +  WioT-i.  '=  \,2,...,SYNPTS 


(4.1) 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


where 

y-i  =  -yi  (4.6) 

and  SIWPTS  is  the  number  of  synthetic  data  points  desired. 
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If  the  test  at  Equation  4.4  fails,  and  if 

lm,ol  <  ,  (4.7) 

then  mirror  imaging  about  the  ocean  surface  is  appropriate,  i.e., 

ciy-i)  =  ciyt),  i  =  IX-.,SYNPTS .  (4.8) 

If  both  of  the  tests  at  Equations  4.4  and  4.7  fail,  then  neither  linear  extrapolation 
nor  mirror  imaging  is  appropriate.  Under  these  circumstances  it  is  appropriate  to  re¬ 
verse  the  rate  of  change  of  the  slope  and  let 

cO_,)  =  c(0)  -f  (m,o  -  /A/«)y_,.,  /  =  \X  -,SYNPTS .  (4.9) 

This  methodology  of  extending  the  real  data  was  implemented  using  five  synthetic 
data  points  both  above  the  surface  and  below  the  bottom.  The  results  are  shown  in 
Tables  20  through  24  and  Figures  28  through  42. 

A  substantial  improvement  in  the  accuracy  of  the  Fourier  series  representation  can 
be  seen  with  respect  to  all  the  functions  examined.  Nevertheless,  some  difficulty  is  still 
evident  near  the  ends  of  the  original  data.  A  final  attempt  to  resolve  this  discrepancy 
was  made  by  using  the  Akima  cubic  spline  representation  of  sound  speed  as  input  to  the 
Fourier  algorithm,  and  sampling  the  SSP  generated  by  the  Akima  cubic  spline  fit  at  some 
arbitrarily  high  frequency.  This  was  attempted  using  a  sampling  rate  ten  times  greater 
than  available  from  the  original  data.  The  synthetic  data  as  described  above  was  not  fit 
using  the  Akima  cubic  spline  method.  Results  for  the  half-period  sine  wave  and 
parabolic  SSPs  are  presented  in  Tables  25  and  26  and  Figures  43  to  48. 

These  results  are  once  again  an  improvement  on  previous  work.  Nonetheless,  some 
inaccuracies  remain.  Unfortunately,  time  constraints  precluded  further,  more  sophisti¬ 
cated  attempts  to  further  improve  the  accuracy  of  the  Fourier  representation. 


60 


Table  20.  FOURIER  SERIES  RESULTS  WITH  LANCZOS 
SMOOTHING  AND  SYNTHETIC  DATA  - 
CONSTANT  SPEED  OF  SOUND 
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LINEAR  SOUND-SPEED  PROFILE  HITIl  A  POSITIVE  GRADIENT 
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Table  22.  FOURIER  SERIES  RESULTS  WITH  LANCZOS 
SMOOTHING  AND  SYNTHETIC  DATA  -  LIN¬ 
EAR  SSP  WITH  NEGATIVE  GRADIENT 
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NUMBER  OF  SYNTHETIC  DATA  POINTS  =  5 

TECHNIQUE:  SPATIAL  FOURIER  SERIES  FIT  HITH  LANCZOS  SMOOTHING 
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SMOOTHING  AND  SYNTHETIC  DATA  - 
HALF-PERIOD  SINE  WAVE  SSP 


»  WCO«HO>A>rO 

* . 

"k  rvKlMOiHCMiH 
♦  O'  .H  III 


. * 

I  I  I  I  I  I  I  I  *-4  A 

I  i  » 

» 


CH  «± 

M 

esi 

M 

M 

<VJ 

Csl 

O 

o 

o 

o 

O 

O 

o 

o 

in 

1 

t 

1 

1 

1 

1 

1 

1 

1 

Ul 

Ul  Ul 

UJ 

Ul 

Ul 

UJ 

Ul 

r 

M  CD 

^  iO  >T 

-o 

CO 

K\ 

s 

to  9^  *& 

ro 

CO 

M 

pH 

>o 

>C 

«0 

(VI  <0 

9 

-• 

Csl 

pH 

CH 

K> 

'O 

'O'  tft 

o 

o 

o 

o 

O 

o 

o 

o 

o 

K) 

M 

(M 

(VI 

CM 

cvi 

(VI 

o 

o 

O 

O 

o 

O 

o 

o 

i/i 

♦ 

1 

1 

1 

i 

1 

1 

1 

Ul 

ur 

Ui 

Ul 

UJ 

Ul 

UJ 

Ul 

r 

o 

o 

>o 

O 

■o 

M 

o 

'O 

N 

o 

iO 

o 

o 

Cl 

'O 

O' 

O' 

o 

H> 

O' 

flO 

>o 

W> 

O' 

— 

o 

O' 

pH 

M 

tA 

O' 

O' 

u> 

o 

o 

o 

o 

O 

O 

o 

o 

ooooooooooooo 

i  I  t  I  I  I  I  I  I  I  I  I  1 

tn  m  n«  ttl  ttt  ti»  ««i  Hi  >«»  in  4I» 

OsO'>OO'O'esfC0t^K)H%vfK^0) 

CnU1^4A(A(A4‘'^mCSi^C^r-l 

ooooooooooooo 


OOOOOOOOOOOOiH 
I  t  I  I  I  I  I  I  i  I  I  I  i 

uj  tn  til  t«t  HI  ttJ  lAi  i«i  tiA  in  >««  HI  >ti 

f^-MO^iAr^'OOCNjNflONOoa) 

>0^>OC^NOO'0'^OCNiOOUIC>J 

00<HOCO'4‘0^tA>OCOO^sOFH 

ooooooooooooo 


SSStUuiS?? 

« 

OOOOOOOO  OpHIAf-^OOOOOOOOO 

oooooooo  ooooooooooooo 

uiUiiLui.iiiiiUjiL  UiUiujuiuiuiuiuiuiuiaiuui 

U)OfOOU)(0<t^MO^>^raO^ 

SSiSSIwJiSjS 

oododdoo  ddddddddddddd 


oooooooo  OMIAMOOOOOOOOO 
oooooooo  OOOOOOOOOOOOO 

Uiuiuiaiuiuiuiui  uiuiujuiuioiuiujujuiuiujuj 
“o'lnSS.^Sm  'tK>wpo^4•Kl«)f^f>.a'lno^^* 
c^S^SSKrsoS  o>to^i-40orvrv^>ors.w 
SoRSor^SS  CV4OOO<MW»AI^O>rK00- 

frtlOMKIMCMCsJO 


OOOOOOOO 
I  I  1  I  I  I  I  I 


OOOOOOOOOOOOO 

I  I  I 


coHio#4inoo<sj  '4'<ntfiu)>^Moou)^<0Hna> 

ooo3SSSf3 


oooooooo 

oooooooo 


ooooooooooooo 

ooooooooooooo 


wR  dddddddd  'eu)uiin^odo<Mdd«M«o^ 
Cr  SSS^SSSSS  f^r^.r'.r-.i^i^oocDOC'O'O' 

y  ^  tm4  ^ 

I-  O^N>P**MMe0r^  ^JrOOM^J^^OK^M^^M^O 

d'ANoinNO'i^  sstsidjsj^sisiisssjss 

lev  O(^o^ol30fla^*•^s 

uj£  s??®  "S55  222i2212£2222 

pH  pH  pH  ^H  pH  pH  pH 


oooooooo  ggggsgggogggg 
oooooooo  ooooooooooooo 

dodddddd  ooooooooooooo 

pHpHpHpHpHpHpHpHpHpHM 


NUHBER  OF  0RI6IKAL  DATA  POINTS  >  11 

NUMBER  OF  SYNTHETIC  DATA  POINTS  <  5 

TECHNIQUE:  SPATIAL  FOURIER  SERIES  FIT  WITH  LANCZ03  SMOOTHING 


Table  24.  FOURIER  SERIES  RESULTS  WITH  LANCZOS 
SMOOTHING  AND  SYNTHETIC  DATA  - 
PARABOLIC  SSP 
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TECHNIQUES  FOR  NUMERICAL  INTERPOLATION  OF  ClY) 


Figure  28.  Sound-Speed  Fit  -  Constant  Speed  of  Sound  with  Synthetic  Data. 
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TECHNIQUES  FOR  NUMERICRL  INTERPOLATION  OF  C(Y) 


Figure  29.  First-Order  Derivative  Fit  -  Constant  Speed  of  Sound  with  Syntlietic 


Data. 
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TECHNIQUES  FOR  NUMERICAL  INTERPOLATION  OF  C(Y) 


Figure  30.  Second-Order  Derivative  Fit  -  Constant  Speed  of  Sound  with  Synthetic 


Data. 
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TECHNIQUES  FOR  NUHERICRL  INTERPOLATION  OF  C(Y) 


Figure  31.  Sound-Speed  Fit  -  Linear  SSP  with  a  Positive  Gradient  and  Synthetic 


Data. 
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TECHNIQUES  FOR  NJHCRICflL  INTERPOLRTION  OF  C(Y) 


Figure  32.  First-Order  Derivative  Fit  -  Linear  SSP  with  a  Positive  Gradient  and 


Synthetic  Data. 
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TECHNIQUES  FOR  NUMERICRL  INTERPOLRTION  OF  C(Y) 


Figure  34.  Sound-Speed  Fit  -  Linear  SSP  «ith  a  Negative  Gradient  and  Synthetic 
Data. 
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TECHNIQUES  FOR  NUHERICRL  INTERPDLRTION  OF  CtY) 


Figure  35.  First-Order  Derivative  Fit  -  Linear  SSP  with  a  Negative  Gradient  and 
Synthetic  Data. 
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TECHNIQUES  FOR  NUMERICRL  INTERPOLRTION  OP  C(Y) 
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TECHNIQUES  FOR  NUMERICRL  INTERPOLnTION  OF  C(Y) 


Figure  38.  First-Order  Derivative  Fit  -  Half-Period  Sine  Wave  SSP  with  Synthetic 
Data. 
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TECHNIQUES  FOR  NUhERICflL  INTCRPDLRTION  OF  C(Y) 


Figure  39.  Second-Order  Derivative  Fit  -  Half-Period  Sine  Wave  SSP  with  Syn¬ 
thetic  Data. 
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TECHNIQUES  FOR  NUHERICRL  INTCRPDLRTION  OF  C(Y) 
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Figure  41.  First-Order  Derivative  Fit  -  Parabolic  SSP  with  Synthetic  Data. 
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TECHNIQUES  FOR  NUMERICAL  INTERPOLATION  OF  CIY) 
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Table  25.  FOURIER  SERIES  RESULTS  WITH  LANCZOS 
SMOOTHING,  SYNTHETIC  DATA,  AKIMA 
CUBIC  SPLINE  INPUT  TO  FOURIER 
METHOD  -  HALF-PERIOD  SINE  WAVE  SSP 
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Table  26.  FOURIER  SERIES  RESULTS  WITH  LANCZOS  SMOOTHING, 
SYNTHETIC  DATA,  AKIM  A  CUBIC  SPLINE  INPUT  TO  FOURIER 
METHOD  -  PARABOLIC  SSP 
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TECHNIQUES  FOR  NUHERlCflL  INTERPOLATION  OF  C(Y) 


Figure  43.  Sound-Speed  Fit  -  Half-Period  Sine  Wave  SSP,  Akinia  Cubic  Spline 


Input  to  Fourier  Method. 
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TECHNIQUES  FOR  NUhERICRL  IKITERPG^RTION  OF  ClY) 


Figure  44.  First-Order  Derivative  Fit  -  Half-Period  Sine  Wave  SSP,  Akima  Cubic 
Spline  Input  to  Fourier  Method. 
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bic  Spline  Input  to  Fourier  Method. 
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Figure  46.  Sound-Speed  Fit  -  Parabolic  SSP,  Akiina  Cubic  Spline  Input  to  Fourier 


Method. 
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TECHNIQUES  FOR  NUMERICAL  INTERPOLATION  OF  CIT) 


Figure  48.  Second-Order  Derivative  Fit  -  Parabolic  SSP,  Akima  Cubic  Spline  Input 
to  Fourier  Method. 
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B.  RAY  TRACE  RESULTS  USING  SPATIAL  FOURIER  SERIES 
REPRESENTATION  OF  SOUND-SPEED  DATA  WITH  THE  RRA  ALGORITHM 
1.  Sound  Speed  Including  Deterministic  Component  Only 

'1  he  RICA  algorithm  using  the  simple  Fourier  series  sound-speed  representation 
including  Lanezos  smoothing  hut  without  synthetic  data  or  Akima  cubic  spline  input  to 
the  Fourier  method  was  tested  using  the  same  source  locations,  SSPs,  and  rays  as  in 
Chapter  III,  Section  A.  That  is.  the  SSPs  were  once  again  defined  as  follows: 

•  SSP  1:  Constant  speed  of  sound  (c  =  ISOOmjs) 

•  SSP  2:  Linear  SSP  with  a  positive  gradient  ( c(^v)  =  c(())  +  gy  m's,  where  the 

speed  of  sound  at  the  surface  c(0)  =  c(>-  =  0)  =  1500  m;s  and  g  =  -f-0.017  I  s  ) 

•  SSP  3;  Linear  SSP  with  a  negative  gradient  ( c(_v)  =  c(0)  -h  gy  m  s.  where  the 

speed  of  sound  at  the  surface  c(0)  =  c(p  =  0)  =  1500  m/s  and  g  =  —0.017  Is)  ■ 

•  SSP  4:  Parabolic  SSP  with  vertex  at  >o  where  cQo)  =  1-^90  m  s  and  the  speed  of 
sound  at  both  the  ocean  surface  and  bottom  is  equal  to  1500  m;s 

(c(j)  =  c(>’3)  +  a(y  -  100)^  m/s  with  a  =  0.001  for  an  ocean  depth  of  200 
m  ). 

The  sound  source  spatial  location  was  established  at  .r,,  =  0  m  and  Zq  =  0  m.  A 
source  depth  ofy^  =  100  m  was  established  for  an  ocean  depth  of  200  m,  unless  otherwise 
noted.  Launch  angles  /?o  w’ere  45'",  85°  and  135“  for  SSPs  1,  2  and  3;  and  S5“,  95“,  and 
135“  for  SSP  4,  which  again  was  designed  to  simulate  a  SOFAR  channel. 

The  ODE  algorithm  using  the  Akima  cubic  spline  fit  to  sound-speed  data  re¬ 
mained  as  the  benchmark  for  comparison  purposes,  as  detailed  in  Chapter  111,  Section 
A.  For  ease  of  reference,  the  results  of  the  ODE  algorithm  are  reproduced  in  the  tables 
of  this  sub-section.  The  ODE's  range  step  size  was  left  at  2.0  m  and  the  arc  length  step 
size  in  the  RRA  algorithm  was  left  at  2.0  m  as  well,  except  where  otherwise  noted. 

The  results  of  the  various  computer  simulation  runs  are  presented  in  1  ables  27 
to  33  below. 
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Table  27.  RAY  ACOUSTICS  RESULTS  FOR  SSP  1,  RAYS  1  THROUGH  3, 
FOURIER  SOUND-SPEED  FIT 


Method 

>0  (m) 

/?Q  (deg) 

r  (km) 

v(m) 

/?  (deg) 

T  (sec) 

s  (km) 

ODE 

100.0 

45.0 

10.0 

100.00 

45.000 

9.428090 

14.1421 

RRA 

100.0 

45.0 

10.0 

99.94 

45.000 

9.42S0S8 

14.1421 

ODE 

100.0 

85.0 

10.0 

174.89 

85.000 

6.692132 

10.0382 

RR.A 

100.0 

85.0 

10.0 

174.88 

85.000 

6.692132 

10.0382 

ODE 

100.0 

135.0 

10.0 

100.00 

135.000 

9.428090 

14.1421 

RRz\ 

100.0 

135.0 

10.0 

100.06 

135.000 

9.428088 

14.1421 

For  SSP  1  (constant  speed  of  sound),  Table  27  shows  that  at  a  horizontal  range 
of  r  =  10  km,  the  differences  between  the  benchmark  ODE  algorithm  and  the  RRA  al¬ 
gorithm  using  the  Fourier  series  sound-speed  fit  are  as  follows: 

•  between  0.01  and  0.06  m  in  depth  and 

•  2  iz  sec  in  travel  time  t. 

Note  the  significant  degradation  of  the  RRA  results  relative  to  those  obtained 
using  the  Akima  cubic  spline  fit  as  in  Table  1  of  Chapter  III,  Section  A.  This  compar¬ 
ative  degradation  of  consistency  between  the  ODE  and  RRA  algorithms  will  be  evident 
throughout  this  subsection. 
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Table  28.  RAY  ACOUSTICS  RESULTS  FOR  SSP  2,  RAYS  I  THROUGH  3, 
FOURIER  SOUND-SPEED  FIT 


Method 

>0  (m) 

/Jo  (deg) 

r  (km) 

y  (m) 

0  (deg) 

T  (sec) 

s  (km) 

ODE 

100.0 

45.0 

10.0 

99.99 

45.000 

9.417413 

14.1421 

RRA 

100.0 

45.0 

10.0 

99.97 

45.000 

9.417429 

14.1421 

ODE 

100.0 

85.0 

10.0 

162.5( 

85.487 

6.683492 

10.0374 

RRA 

100.0 

85.0 

10.0 

163.31 

85.482 

6.683550 

10.0374 

ODE 

100.0 

135.0 

10.0 

100.01 

135.000 

9.417413 

14.1421 

RRA 

100.0 

85.0 

10.0 

100.03 

135.000 

9.417430 

14.1422 

For  SSP  2  (linear  SSP  with  a  positive  gradient),  Table  28  shows  that  at  a  hori¬ 
zontal  range  r  =  10  km,  the  diflerences  between  the  benchmark  ODE  algorithm  and  the 
RRA  algorithm  using  the  Fourier  series  sound-speed  fit  are  as  follows: 

•  between  0.02  and  0. 19  m  in  depth  y, 

•  0.005  deg  in  angle  of  arrival 

•  between  16  and  58  n  sec  in  travel  time  t,  and 

•  0. 1  m  in  path  length  s. 


I 
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Table  29.  RAY  ACOUSTICS  RESULTS  FOR  SSP  3,  RAYS  1  THROUGH  3, 
FOURIER  SOUND-SPEED  FIT 


Method 

>0  (m) 

/^O  (deg) 

r  (km) 

y  (m) 

/?  (deg) 

T  (sec) 

s  (km) 

ODE 

100.0 

45.0 

10.0 

99.99 

45.000 

9.438784 

14.1421 

RRA 

100.0 

45.0 

10.0 

99.96 

45.000 

9.438799 

14.1422 

ODE 

lOO.O 

85.0 

lO.O 

168.89 

84.512 

6.699404 

10.0380 

RRA 

100.0 

85.0 

— 

10.0 

169.76 

84.492 

6.699471 

10.0380 

ODE 

100.0 

135.0 

10.0 

100.01 

135.000 

9.438784 

14.1421 

RRA 

100.0 

85.0 

10.0 

100.04 

135.000 

9.438798 

14.1421 

For  SSP  3  (linear  SSP  with  a  negative  gradient),  Table  29  shows  that  at  a  hori¬ 
zontal  range  10  km,  the  dilTerence  between  the  benchmark  ODE  algorithm  and  the 
RRA  algorithm  using  the  Fourier  sound-speed  fit  are  as  follows: 

•  between  0.03  and  0.13  m  in  final  depth 

•  0.020  deg  in  angle  of  arrival  /?, 

•  between  14  and  67  n  sec  in  travel  time  t,  and 

•  0. 1  m  in  path  length  s. 

The  results  contained  in  the  above  tables,  while  poorer  than  those  detailed  in 
Chapter  III,  Section  A,  nevertheless  produce  travel  times  which  are  consistent  in  the 
millisecond  region.  However,  this  is  not  the  case  for  the  SSPs  e.xamined  next. 
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Table  30.  RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3, 
FOURIER  SOUND-SPEED  FIT 


Method 

To  (m) 

/Jo  (deg) 

r  (km) 

T  (m) 

/?  (deg) 

r  (see) 

5  (km) 

ODE 

100.0 

45.0 

10.0 

37.06 

87.244 

6.710387 

10.0183 

RRA 

lOO.O 

45.0 

10.0 

42.86 

86.674 

6.709733 

10.0185 

ODE 

100.0 

85.0 

10.0 

162.94 

92.756 

6.710387 

10.0183 

RR.\ 

100.0 

85.0 

10.0 

157.14. 

93.326 

6.709733 

10.0185 

ODE 

100.0 

135.0 

10.0 

144.97 

134.922 

9.448852 

14.1104 

RR.\ 

100.0 

85.0 

10.0 

143.98 

134.930 

9.448515 

14.1112 

For  SSP  4  (parabolic  SSP),  Table  30  shows  that  at  a  horizontal  range  of  r  =  10 
km,  the  differences  between  the  benchmark  ODE  algorithm  and  the  RRA  algorithm 
using  the  Fourier  series  sound-speed  fit  are  as  follows: 

•  between  0.99  and  5.8  m  in  final  depth 

•  between  0.008  and  0.570  deg  in  angle  of  arrival  /?, 

•  between  337  and  654  n  sec  in  travel  time  t,  and 

•  0.2  m  in  path  length  s. 

As  in  Chapter  111,  Section  A,  SSP  4  shows  significantly  less  consistency  be¬ 
tween  the  two  algorithms  than  the  previous  SSPs  examined.  Again,  two  approaches 
were  taken  to  attempt  to  reduce  this  discrepancy.  First,  the  arc  length  step  size  was  re¬ 
duced  from  2.0  to  1.0  m  (see  Tables  12  and  13)  and  second,  the  steepness  of  the  gradient 
was  reduced  from  0.1  to  0.02  l/sec  (see  Table  14). 


Table  31.  RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  2.0 
M  ARC  LENGTH  STEP  SIZE,  FOURIER  SOUND-SPEED  FIT 


Method 

>0  (m) 

/?0  (deg) 

r  (km) 

y  (ni) 

/?  (deg) 

T  (sec) 

s  (km) 

ODE 

100.0 

45.0 

5.0 

64.25 

S5.591 

3.354779 

5.0088 

RR.\ 

100.0 

45.0 

5.0 

68.79 

85.435 

3.354591 

5.0090 

ODE 

100.0 

85.0 

5.0 

135.75 

94.403 

3.354779 

5.0088 

RRA 

100.0 

85  0 

5.0 

131.21 

94.565 

3.354591 

5.0090 

ODE 

100.0 

135.0 

5.0 

77.50 

45.019 

4.724411 

7.0552 

RR.\ 

100.0 

85.0 

5.0 

77.82 

45.017 

4.724162 

7.0554 

For  SSP  4  (parabolic  sound  speed),  Table  31  shows  that  at  a  horizontal  range 
of  r  =  5  km,  using  an  arc  length  step  size  of  2.0  m,  the  differences  between  the 
benchmark  ODE  algorithm  and  the  RRA  algorithm  using  the  Fourier  series  sound-speed 
fit  are  as  follows; 

•  between  0.32  and  4.54  m  in  final  depth 

•  between  0.002  and  0.162  deg  in  angle  of  arrival  /?, 

•  between  188  and  249  n  sec  in  travel  time  t,  and 

•  0.2  m  in  path  length  s. 

The  data  in  Table  31  above  was  generated  in  order  to  enable  a  direct  compar¬ 
ison  with  the  arc  length  step  size  set  to  1.0  m,  as  shown  in  Table  32  below.  The  re¬ 
duction  in  final  range  from  10  to  5  km  was  necessitated  by  the  structure  of  the  RRA 
computer  code,  which  required  large  quantities  of  computer  storage  space  to  be  allo¬ 
cated,  particularly  as  arc  length  step  size  was  reduced. 


Table  32.  RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  1  THROUGH  3,  1.0 
M  ARC  LENGTH  STEP  SIZE 


Method 

>0  (ni) 

po  (deg) 

— 

r  (km) 

— 

y  (m) 

P  (deg) 

T  (sec) 

s  (km) 

ODE 

100.0 

45.0 

5.0 

64.25 

85.597 

3.354779 

5.0088 

RRA 

100.0 

45.0 

5.0 

68.86 

85.442 

3.354604 

5.0090 

ODE 

100.0 

85.0 

5.0 

135.75 

94.403 

3.354779 

5.0088 

RRA 

100.0 

85.0 

5.0 

131.14 

94.558 

3.354604 

5.0090 

ODE 

100.0 

135.0 

5.0 

77.50 

45.019 

4.724411 

7.0552 

RRA 

100.0 

85.0 

5.0 

ll.lA 

45.019 

4.724126 

7.0554 

For  SSP  4  (parabolic  SSP),  Table  32  shows  that  at  a  horizontal  range  orr=  5 
km,  and  with  the  RRA  arc  length  step  size  set  to  1.0  m,  the  difierences  between  the 
benchmark  ODE  algorithm  and  the  RRA  algorithm  using  the  Fourier  series  sound-speed 
fit  are  as  follows; 

•  between  0.24  and  4.61  m  in  final  depth 

•  0.155  deg  in  angle  of  arrival /J, 

•  between  175  and  285  fi  sec  in  travel  time  r,  and 

•  0.2  m  in  path  length  s. 
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Table  33.  RAY  ACOUSTICS  RESULTS  FOR  SSP  4,  RAYS  I  THROUGH  3,  2.0 
M  ARC  LENGTH  STEP  SIZE,  BOTTOM  DEPTH  1000  M,  FOURIER 
SOUND  SPEED  FIT 


Method 

>0  (m) 

Ao (deg) 

r  (km) 

y  (ni) 

P  (deg) 

T  (sec) 

s  (km) 

ODE 

100.0 

45.0 

10.0 

777.41 

93.388 

6.705898 

10.0150 

RRA 

100.0 

45.0 

10.0 

768.46 

93.522 

6.705420 

10.0152 

ODE 

100.0 

85.0 

10.0 

222.59 

86.612 

6.705898 

10.0150 

RRA 

100.0 

85.0 

10.0 

231.54 

86.478 

6.705420 

10.0152 

ODE 

100.0 

135.0 

10.0 

545.01 

134.997 

9.448829 

14.1104 

RRA 

100.0 

85.0 

10.0 

544.58 

134.997 

9.448180 

14.1107 

For  SSP  4  (parabolic  speed  of  sound),  with  vertex  c(5(X))  =  1490  m,'s  and  speed 
of  sound  at  the  ocean  surface  c(0)  =  1500  m/s  and  speed  of  sound  at  the  ocean  bottom 
c(lOOO)  =  1500  m,'s,  i.e.,  an  ocean  depth  of  1000  m;  at  a  horizontal  range  of  r  =  10  km, 
the  diflerences  between  the  benchmark  ODE  algorithm  and  the  RRA  algorithm  using 
the  Fourier  series  sound-speed  fit  are  as  foUows: 

•  between  0.43  and  8.95  m  in  final  depth  y, 

•  0.134  deg  in  angle  of  arrival 

•  between  478  and  649  n  sec  in  travel  time  t,  and 

•  between  0.2  and  0.3  m  in  path  length  s. 

Table  33  shows  no  significant  improvement  over  Table  30  for  SSP  4  with  the 
original  extreme  gradient.  The  errors  associated  with  the  Fourier  sound  speed  repre¬ 
sentation  have  resulted  in  this  lack  of  improvement,  whereas  substantial  improvement 
was  observed  between  these  two  cases  for  the  Akima  cubic  spline  fit  with  the  RRA  al¬ 
gorithm  (see  Tables  4  and  7  in  Chapter  HI,  Section  A). 
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2.  Sound  Speed  Including  Deterministic  and  Random  Components 

In  order  to  investigate  the  efTects  of  randomness  in  the  SSP,  a  mod(^\  of  the 
power  spectral  density  of  the  random  component  of  the  speed  of  sound  had  to  be 
proposed.  The  most  straightforward  model  of  band-limited  white  noise  was  used,  as 
given  by 


l/vl  <  SIV 

141  >  BIV 


(4.10) 


where  BW  is  the  bandwidth,  and  5^  is  the  constant  magnitude  of  the  power  spectrum  at 
all  spatial  frequencies  4  • 

As  a  consequence  of  the  assumption  of  zero  mean  white  noise,  the  variance  of 
is  given  by 

<7^^  =  2SoBfV  (4.11) 

and,  as  a  result, 

A  value  for  the  variance  of  c^ly)  in  shallow  water  of  2.5  x  10-’  was  obtained 
from  Brekhovskikh  and  Lysanov  [Ref  8:  p.  208  |.  The  simulation  program  was  run  with 
randomness  incorporated  as  above  and  for  SSP  1  (constant  speed  of  sound).  When  the 
results  were  compared  with  those  from  the  RRA  algorithm  previously  run  using  the 
Fourier  series  representation  of  sound-speed  data,  differences  of  a  maximum  of  0.1  m  in 
final  depth  y,  0.002  deg  in  angle  of  arrival  ft,  and  2  fi  sec  in  travel  time  t  were  observed. 
The  RRA  algorithm  was  then  run  with  SSP  4  (parabolic  SSP  in  a  200  m  deep  ocean) 
incorporating  randomness  and  differences  of  a  maximum  of  0.22  m  in  final  depth  y  , 
0.018  deg  in  angle  of  arrival  ft,  and  16  /x  sec  in  travel  time  t  were  observed. 

The  variance  was  then  increased  to  determine  at  what  magnitude  it  caused  the 
randomness  to  become  significant.  Further  simulation  runs  revealed  that  for  SSP  1, 
differences  in  travel  times  of  greater  than  10  n  sec  needed  al^  to  be  at  least  of  the  order 
of  10~’  ,  four  orders  of  magnitude  greater  than  the  estimate  given  in  Brekhovskikh  and 
Lysanov  [Ref  8:  p.  208  ]. 


2BIV 


(4.10) 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  recursive  ray  acoustics  (RRA)  algorithm  was  showrt  to  be  an  accurate  ray 
acoustics  algorithm  when  used  in  conjunction  with  an  Akima  cubic  spline  fit  to  sound- 
speed  data.  The  RRA  algorithm  can  be  made  as  accurate  as  desired  by  reducing  the  arc 
length  step  size.  Although  the  ODE  algorithm  computes  amplitudes  which  the  RR,\ 
algorithm  does  not,  the  RRA  algorithm  was  observed  to  be  significantly  faster  in  terms 
of  computer  execution  time. 

The  ability  of  the  RR.A  algorithm  to  handle  sound  speed  as  a  function  of  all  three 
spatial  coordinates  enables  its  application  to  real  world  problems  wherein  a  horizontal 
gradient  of  sound  speed  exists,  either  in  the  cross-range  or  down-range  direction.  This 
could  conceivably  be  the  case  for  either  long  range  ray  propagation  or  for  propagation 
in  estuaries  or  across  current  boundaries  where  significant  horizontal  gradients  can  exist. 

A  spatial  Fourier  series  (SFS)  representation  of  sound-speed  data,  while  attractive 
for  reasons  of  its  spectral  content  and  the  ease  of  incorporation  of  randomness  at  spe¬ 
cific  harmonics,  encountered  dilTiculty  with  regard  to  providing  accurate  fits  to  sound- 
speed  data.  Several  techniques  were  used  to  attempt  to  rectify  this,  including  the  use 
of  Lanezos  smoothing  ,  addition  of  synthetic  data  above  the  ocean  surface  and  below 
the  ocean  bottom,  and  sampling  of  the  Akima  cubic  spline  fit  to  provide  data  points  for 
input  to  the  SFS  method.  These  techniques  resulted  in  some  substantial  success  in  the 
reduction  of  the  inaccuracies  of  the  SFS  method.  However,  although  the  SE'S  method 
eventually  provided  accurate  fits  to  the  sound-speed  data,  it  still  provided  inaccurate  first 
and  second-order  derivatives  when  compared  against  the  results  obtained  from  the 
Akima  cubic  spline  method. 

A  comparison  was  made  between  the  ODE  algorithm  using  an  Akima  cubic  spline 
sound-speed  representation  and  the  RRA  algorithm  using  a  SFS  sound-speed  represen¬ 
tation,  with  the  SFS  representation  in  its  most  basic  form,  that  is,  with  none  of  the 
techniques  referred  to  in  the  above  paragraph  incorporated.  The  RRA/SFS  combina¬ 
tion  was  found  to  be  w'anting  in  accuracy. 

As  was  mentioned  previously,  the  RRA  algorithm  using  the  Akima  cubic  spline  fit 
to  sound-speed  data  was  proven  fast  and  accurate,  and  can  be  made  as  accurate  as  de¬ 
sired  by  reducing  the  arc  length  step  size.  However,  investigation  is  warranted  to  de- 
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termine  if  the  RRA  algorithm  can  be  generalized  to  include  amplitude  calculations, 
especially  at  focal  points  and,  hence,  caustics. 

Throughout  the  simulations  undertaken  in  this  thesis,  a  zero-slope,  flat  ocean  sur¬ 
face  and  bottom  was  assumed.  In  order  to  apply  the  RRA  algorithm  to  realistic  ocean 
boundaries,  the  ability  to  accommodate  arbitrarily  shaped  surfaces  should  be  included. 
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